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ABSTRACT 


A  technique  Is  presented  for  the  constructions  of  stress* 
strain  relations  based  on  experimental,  cyclic,  damping  data. 

The  extension  of  this  technique  to  the  biaxial  stress  case  is 
then  shown  followed  by  an  exai^le  application  involving  flexural 
vibrations  of  a  cantilever. 
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I.  INTRODUCTION 


The  purpose  of  this  paper  Is  to  develop  and  demonstrate  a 
simple  mathematical  formulation  of  one  and  two  dimensional,  damped 
stress-strain  relations  based  on  a  special  reduction  of  experi¬ 
mental  data.  The  motivation  for  this  work,  which  might  be  consid¬ 
ered  as  an  Intermediate  approach,  stems  from  the  marked  complexity 
and  large  number  (at  least  a  dozen)  of  separate  physical  phenomena 
which  have  been  discovered  as  being  integrant  to  material  damp¬ 
ing.  (1)*  The  implication  is  that  a  rigorous  approach,  which  per¬ 
force  must  reflect  svich  involvement,  could  only  be  speculative 
with  the  present  limited  experimental  data.  In  spite  of  the 
physical  complexity  of  material  damping,  however,  certain  system¬ 
atizations  of  the  experimental  data  have  been  achieved,  and  it  is 
for  such  a  development  that  the  present  paper  constructs  a  mathe¬ 
matical  theory. 

The  simplest  mathematical  formulation  for  energy  dissipation 
which  can  be  used  to  describe  structviral  damping  is  that  which 
corresponds  to  the  linear  dashpot.  Its  plot  on  stress-strain 
coordinates,  for  sinusoidal  straining  action,  produces  smooth 
contoxired  hysteresis  loops  as  shown  in  Figure  la.  If  this  formu¬ 
lation  is  moiified  by  dividing  the  damping  coefficient  by  the 
frequency  (implying  restriction  of  attention  to  discrete  fre¬ 
quencies),  then  this  removes  the  linear  frequency  dependence  of 
the  energy  dissipation  per  cycle.  This  is  usually  termed  "hys- 
teretic  damping,"  and  is  one  of  the  most  prominent  characteristics 
of  material  damping.  Thus  for  the  uniaxial  case,  the  stress- 
strain  relation  for  hysteretic  damping  becomes 


^  (1-1) 

where  <t  and  e  are  respectively,  the  stress  and  strain,  E  is 
Young's  modulus,  C  is  the  damping  coefficient,  w  is  the  fre¬ 
quency  of  the  sinusoidal  straining  action  and  i  is  the  strain 
rate. 

Although  it  turns  out  that  the  foregoing  formulation,  which 
produces  elliptical  hysteresis  loops,  can  produce  excellent  re¬ 
sults  in  most  vibration  analysis  problems,  the  experimentally 
obtained  loops  for  metals  have  sharp  comers  at  both  ends,  and 


*Number  in  parentheses  refer  to  a  bibliography  on  page  26. 
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furthermore,  are  practically  linear  in  the  central,  or  low  stress 
strain  region.  This  suggests  that  the  next  stage  of  complexity 
in  the  analytical  development  might  be  an  attempt  at  a  direct, 
geometric  description  of  the  hysteresis  loops  incorporating  the 
aforementioned  characteristics.  This  has  been  done  in  various 
ways.  A  formulation,  due  to  Davidenkov  (2),  which  brings  in  the 
characteristic  sharp  comers  as  shown  in  Figure  lb,  and  which 
also  incorporates  certain  dependence  upon  strain  amplitude,  is 
represented  by  the  following  equations: 


The  two  equations,  one  each  for  the  upper  and  lower  curves  of  the 
hysteresis  loop,  are  now  a  necessary  complication.  The  constants 
/X  and  m  are  left  to  be  evaluated  experimentally  for  the 
material  under  consideration,  and  is  the  strain  amplitude. 

This  particular  description  of  material  damping  has  been  applied 
to  a  number  of  problems  by  Pisarenko  (3) .  The  associated  mathe¬ 
matical  problem,  of  course,  is  nonlinear,  but  his  results  show 
that  a  first  approximation  only,  derived  by  the  use  of  asymptotic 
expansions,  gives  a  fairly  accurate  result. 

Another  example  of  the  foregoing  technique  of  geometric 
description  of  hysteresis  loops  is  that  developed  by  Rang  (4) . 

In  this  case  a  cubic  polynomial  is  used  to  describe  each  of  the 
two  sides  of  the  loop,  which  produces  sharp  comers,  with  the 
constants  being  determined  by  assigning  values  to  certain  charac¬ 
teristic  tangents  of  the  loop  and  the  locus  of  the  comer  points. 
Both  of  the  foregoing  formulations  allow  close  approximation  of 
experimentally  observed  physical  phenomena,  but  at  the  expense  of 
making  subsequent  mathematical  analysis  quite  difficult.  For  the 
Rang  case,  an  example  solution  is  given  only  for  a  discrete,  one 
degree  of  freedom,  spring -mass  model.  Furthermore,  both  of  the 
foregoing  analyses  are  developed  strictly  for  the  uniaxial  stress 
case  and  no  attempt  appears  to  have  been  made  at  incorporating 
them,  in  any  systematic  way,  with  the  biaxial  stress  case  or  with 
other  phenomena  such  as  hereditary  effects,  etc.  All  descriptions, 
however,  can  be  modified  to  include  various  types  of  dependence 
on  strain  amplitude. 

Thus,  at  present,  there  is  not  only  the  penalty  of  mathemat¬ 
ical  complication  in  employing  relatively  precise  descriptions  of 
the  hysteresis  loop,  but  the  added  disadvantage  that  no  suitable 
collection  of  data  is  available  for  a  wide  variety  of  materials. 
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In  fact,  the  only  general  technique  which  has  been  shown  up  to 
now,  for  systematizing  such  data,  uses  the  area  of  the  hysteresis 
loop,  plus  certain  fatigue  properties,  with  no  reference  to  loop 
shape.  In  other  words,  it  does  not  contain  information  with  re¬ 
spect  to  instantaneous  rates  of  energy  dissipation,  which  would 
therefore  have  to  be  added  in  some  way  if  the  more  precise  de¬ 
scriptions  were  to  be  applied.  Besides  this  problem  of  expanding 
the  collected  data,  there  occurs  the  still  additional  problem  of 
interpreting  equivalent  damping  conditions  under  uniaxial  and 
biaxial  stress  conditions.  The  precise  phenomena  involved  in  this 
connection  are  not  fully  understood  at  the  present  time.  The  only 
relevant  theory  for  associating  data  obtained  using  uniaxial 
stress  to  the  biaxial  stress  case  is  that  due  to  Robertson  and 
Yorgiadis  (5),  who  define  equivalent  states  (of  material  damping) 
as  those  which  produce  the  same  distortion  energy,  and  that  due 
to  Mentel  (6),  who  suggests  that  in  such  cases  as  plate  vibration, 
some  fractional  part  of  the  dilatational  strain  energy  ought  to 
be  included  in  the  foregoing  equivalence.  The  latter  formulation 
can  be  expressed  analytically 


="0  =  J  3  T-m  +  A  <7/  (1.3) 

where  is  the  maximiim  (normal)  stress  in  the  uniaxial  case 

and  Tpp  and  are  the  maximiim  shear  and  the  "mean  normal"  or 

dilatational  components  of  stress  in  the  equivalent  biaxial  case. 
The  development  of  this  equation  is  worth  special  notice.  Con¬ 
sider  the  two  states  of  stress  as  shown  in  Figure  2.  If  we 
observe  that  the  principal  stresses  in  the  biaxial  stress  state 
can  be  designated  by 


then  the  expressions  for  the  distortional  and  dilatational  strain 
energies  take  on  the  following  form: 

for  the  uniaxial  stress  state. 
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for  the  biaxial  stress  state, 
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Thus,  if  we  choose  to  ignore  the  mean  normal  stress  ,  say 

by  setting  cr^  =  0,  and  then  equate  the  distort ional  strain 
energies  in  tne  two  cases,  we  obtain 


+  1/  2 
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which  is  the  Robertson-Yorgiadis  criterion.  On  the  other  hand, 
if  we  retain  all  quantities  and  equate  the  total  strain  energies 
in  the  two  cases,  we  find  the  equivalent  expression  to  be 


a  =  ^  2  (  I  +  1/ )  T  *  +  2  (  I  —  1/ )  (T 
n  V  m  d 

If  we  now  set  v  “1/3,  then  (1-5)  becomes 


(1-5) 
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Alternatively,  if  regard  the  damping  phenomenon  as  being  associated 
primarily  with  a  "plastic  type'  straining  process  wherein 
we  obtain 


3t^  + 
m 


(1-7) 


The  coefficient  of  given  by  (1-6)  and  (1-7),  is  thus 

observed  to  be  at  roost  only  moderately  altered  from  its  value  in 
(1-4).  Hence,  in  order  to  allow  for  a  variable  damping  influence 
between  distort ional  and  dilatational  straining,  the  form  given 
by  (1-3)  is  suggested  as  a  simple,  one  parameter  description  for 
empirical  reduction  of  experimental  data.  Although  the  uniaxial 
to  biaxial  damping  equivalence  given  by  (1-3)  carries  the  impli¬ 
cation  that  the  shear  stress  criterion,  as  expressed  by  (1-4), 
gives  the  lower  boxind  for  material  damping,  it  should  be  observed 
that  no  argument  is  given  to  exclude  negative  values  of  A  which 
means,  for  example,  that  a  nullifying  type  of  interaction  could 
be  described.  Note  that  the  first  stress  invariant  I^,  for  the 
biaxial  case,  is  given  by 


I.  =  I 


ij 


=  2  0-. 


so  that  (1-3)  can  be  rewritten 


3r"  ^AI. 


m 


(1-8) 


This  reduction  can  be  carried  one  step  further  by  introducing 
the  second  stress  invariant 


= 


= 


m 


which  gives 


m 


=  —  I*  -  I 

4  I  -^2 


so  that  (1-8)  can  be  further  modified  to 


=  y -^(3+A)If  -  3  1, 


(1-9) 
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It  is  observed,  however,  that  I2  introduces  no  new  functional 
forms  beyond  those  contained  in  ^1-8^,  so  that  a  simple,  one 
parameter  description,  based  on  (1-8),  is 


=  y  ^  A*  I 


(1-10) 


which  is  the  expression  used  in  this  paper.  A  two  parameter 
description,  based  on  (1-9),  would  be 


(1-11) 


This  latter  formulation  might  be  the  most  promising  for  the  tri- 
axial  stress  case  but  this  problem  is  left  for  future  studies. 


In  view  of  the  foregoing  unclear  situation  concerning  mate¬ 
rial  damping  under  biaxial  stress,  it  was  concluded  that  a  useful 
interim  step  for  theoretical  response  work,  would  be  to  construct 
a  mathematical  theory  which  was  counterpart  only  to  the  degree  of 
systematization  thus  far  shown  for  experimental  data.  This  implies 
neglect  of  loop  shape,  so  that  the  simplest  mathematical  shape, 
viz.,  the  ellipse,  might  as  well  be  taken.  The  chief  effect  of 
this  simplification,  whose  primary  geometric  characteristic  is  the 
elimination  of  the  sharp  corners  at  the  loop  tips,  might  be  ex¬ 
pected  to  be  significant  only  at  the  highest  modes  and  frequencies 
of  vibration.  No  such  effect,  however,  is  demonstrated  since  the 
analysis  is  carried  out  only  for  vibrations  in  the  neighborhood  of 
the  fundamental  resonance.  In  this  case,  when  the  loops  corre¬ 
spond  to  structural  metals  at  stress  amplitudes  below  yield,  the 
area  only  and  not  the  loop  shape  appears  to  be  significant  [3] . 
Further  discussion  of  this  question  is  included  in  subsequent 
sections  of  this  paper. 

The  one  parameter  equivalence  relation  between  uniaxial  and 
biaxial  stress  states,  as  expressed  by  (1-10)  and  the  experimental 
uniaxial  damping  data  are  then  used  as  a  basis  for  generating  the 
stress  strain  relations  for  the  biaxial  case.  These  relations 
permit  the  formulation  of  the  vibration  response  problem  for 
plates  exhibiting  material  damping.  Example  solutions  for  the  case 
of  a  circular  plate  are  given  in  [7] . 
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II.  ANALYTICAL  FORMULATION  OF  MATERIAL  DAMPING 


2.1  General  Problem 

It  was  observed  in  the  foregoing  section,  that  at  least  a 
dozen  individual  phenomena  (e.g.,  magnetostriction,  intercrystal 
thermal  currents,  interstitial  solute  atoms,  etc.)  have  been 
identified  as  being  significant  in  the  production  of  material 
damping.  Any  widely  applicable  reduction  of  danroing  data  can 
therefore  be  expected  to  be  either  extremely  elaborate,  which 
singles  out  the  aforementioned  phenomena  and  their  associated 
parameters,  or  extremely  simple,  which  suppresses  them.  At  the 
present  time,  only  the  latter  approach  has  led  to  significant 
success  based  on  a  broad  assimilation  of  data.  This  development 
has  been  due  to  Lazan  who  has  shown,  empirically,  that  the 
fatigue  strength  of  material  can  be  used  as  an  important  parameter 
in  the  reduction  of  experimental  damping  data  (see  Figure  2.10 
of  Ref.  (8)).  This  parameter  is  the  ratio  of  the  cyclic  stress 
Sq  (in  mathematical  notation,  the  stress  amplitude  )  to  the 
fatigue  Strengths^,  corresponding  to  2  x  107  cycles,  both  quan¬ 
tities  corresponding  to  uniaxial,  sinusoidal,  constant  amplitude 
and  fixed  frequency  conditions.  The  experimental  results  show, 
that  for  the  given  selection  of  structural  materials,  the  energy 
dissipation  is  essentially  independent  of  both  the  frequency  of 
cycling  £ind  the  stress  history  when  the  stress  amplitude  is  kept 
below  a  certain  critical  value.  This  observation  confirms  the 
common  notion  of  hysteretic  damping,  which  serves  adequately  in 
many  simplified  analysis  of  vibration  response. 

The  aforementioned  critical  value  of  stress  amplitude  is 
called  the  cyclic  stress  sensitivity  limit  S|_ .  Both  the  fre¬ 
quency  of  cycling  and  stress  history  develop  marked  influence  for 
cyclic  stressing  above  this  stress  limit.  An  empirical  determi¬ 
nation  of  this  critical  value,  due  to  Lazan  (1),  is  80  percent  of 
the  fatigue  strength  corresponding  to  2  x  107  cycles.  Fortunately, 
from  the  analytical  point  of  view,  this  value  is  relatively  high 
compared  to  stress  levels  nominally  encountered  in  structures 
subjected  to  long  term  dynamic  loading. 

The  conclusion  derived  from  the  foregoing  sections  is  that 
at  stress  amplitudes  below  the  cyclic  stress  sensitivity  limit, 
the  cyclic  material  damping  is  adequately  expressed  merely  as  a 
function  of  the  strer  *,  amplitude  Sq.  Such  a  description  has  been 
found  by  Lazan  (8)  to  be 


Do  = 


(2.1-1) 
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where  J  and  a  are  adjusted  for  the  given  material,  and  D-  is  the 
cyclic  energy  dissipation  per  unit  volxaoe.  The  range  of  values 
taken  by  the  parameters  J  and  of  is  relatively  small  (no  changes 
of  order  of  magnitude)  for  a  large  selection  of  structural  mate¬ 
rials.  The  total  cyclic  energy  dissipation  D^,  in  a  given  volume 
of  material  is  then  given  by 


D 


m 


=  / 


Volume 


(2.1-2) 


For  reasons  which  will  be  apparent  later,  the  alteimative  expres 
sion 


D 


0 


=  I 

n 


Jn  s! 


(2.1-3) 


where  n  is  restricted  to  being  an  integer,  will  be  introduced  as 
a  replacement  of  (2.1-1). 


2.2  Uniaxial  Stress 


We  have  already  introduced  the  stress  strain  relations  of 
Davidenkov  as  equations  (1-2)  in  section  1.  If  these  equations 
are  integrated  over  a  complete  cycle  in  order  to  obtain  the  net 
loop  area,  and  hence  the  cyclic  energy  dissipation  Dri*  we  find 
that 


^0  /• "  ^0 

D  =J  o^c/c+J  <T  d  e  = 


.m  +  i 


2""'  (  m  -\)fiE  m+i 


We  now  note  that 


S.  =  =  E 


u.  m- 1 

(  2  €„  ) 

m  0 


}. 


so  that  if  we  assume  small  strains  (and  small  nonlinearity)  in 
the  sense  that 

fji  ,  m-i 

«  I 

we  can  rewrite  equation  (2.2-1)  in  the  form 


2^*'  (  m  -  I  )  ^  ' 

m(m+l)  e"'  ° 


(2.2-2) 
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Equation  (2.2-2)  is  now  in  the  same  form  as  the  empirical  expres¬ 
sion  (2.1-1)  and  we  obtain  an  identity  if  we  choose 

m  =  a  —  I 


and 


H-  = 


Q>(o!-|)E°^  ' 
2°^  {ci-2) 


(2.2-3) 


The  fact  of  being  able  to  demonstrate  the  foregoing  equivalence 
lends  additional  validity  both  to  the  experimentally  derived 
equation  (2.1-1)  and  the  Davidenkov  equations  (1-2).  However, 
the  usefulness  of  the  Davidenkov  equations  is  limited  by  the 
complication  of  having  to  keep  track  of  the  proper  branch  of  the 
hysteresis  loop  for  different  material  points  in  addition  to 
having  mathematical  nonlinearity.  Furthermore,  the  value  of  OL 
in  (2.1-1)  lies  between  2  and  3  for  most  structural  metals  with 
the  value  being  nearer  2  for  most  steels  (8  )  (  a  *=  2  for  1020 
steel,  ct  =2.3  for  Sandvik  steel  Q-T) .  This  corresponds  to 
m  =  1  in  (1-2),  for  which  value  these  equations  collapse  to 


a  =  E  {  I  —  /i.  )  € 

which  describes  simply  the  undamped  elastic  case.  Thus,  for 
integral  values  of  d  ,  equations  (1-2)  are  useful  only  for  the 
extreme  case  of  a  =  3  if  material  damping  is  to  be  included. 

The  approach  to  the  development  of  damped  stress-strain  rela¬ 
tions  proposed  in  this  paper,  is  to  use  strain  rate  as  the 
device  to  introduce  elliptical  hysteresis  loops  (in  the  manner  of 
complex  damping  [9,  10]  ).  The  advantage  of  this  approach  is  that 
it  allows  direct  and  simple  usage  of  the  full  hysteresis  loop. 

The  stress-strain  relation  can  thus  be  written  in  the  form 


where  Ci  and  /3  account  for  elastic  nonlinearity  and  C2  accounts 
for  the  damping.  The  available  experimental  data,  however,  deals 
primeirily  with  structural  metals,  where,  in  the  range  of  stress 
levels  at  which  (2.1-1)  is  applicable,  the  elastic  component  is 
strictly  linear.  Hence,  in  place  of  (2.2-4),  we  are  able  to  use 
the  simpler  relation 


<J  =  E  {  €  +  C  €  ) 


(2.2-5) 


9 


where  the  damping  coefficient  C  is  expected  to  be  a  function  of 
the  strain  amplitude,  frequency  and  the  elastic  constants. 

Integrating  (2,2-5)  for  the  cyclic  damping,  we  obtain 


D.  =  E  C 


/, 


Cycle 


(2.2-6) 


If  the  straining  action  is  now  assumed  to  have  the  sinusoidal 
variat ion 


€  =  cos  t  {  >  0  ) 

the  stress  strain  relation  gives 

c  =  E  6o  ^  cos  ojg  t  —  C  sin  tUg  t  )  , 


so  that 


So  =  E  (  I  +  C""  Wg  )  2 

the  energy  dissipation,  as  expressed  by  (2.2-6), 

27r 

Do  =  =  TTECcUg  6^" 

substitution  of  (2.2-7)  into  (2.2-1)  gives 


(2.2-7) 
then  becomes 


(2.2-8) 


^0  ^  e;  Cq  (  I  +  c  o^g ) 


(2.2-9) 


which  can  be  compared  with  (2.2-8).  Requiring  the  equivalence  of 
(2.2-8)  with  (2.2-9)  thus  gives 


TT  Co^g  =  J  E®  '  *  (  I  +  CUg  )  ® 


(2.2-10) 


Equation  (2.2-10)  identifies  the  functional  dependence  required  of 
the  damping  parameter  C  •  This  expression  can  be  simplified  if 
we  observe  that,  for  most  structural  applications,  we  can  expect 
<<  I  .  The  resulting  expression  for  C  is 
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c 


(2.2-11) 


J  p®”' 

TT^s  ^ 

Recalling  that  for  structural  metals  a  usually  has  some  non¬ 
integral  value  between  2  and  3,  the  alternative  expression  (2.1-3) 
can  be  employed  to  retain  the  simplicity  of  being  able  to  stay 
with  integers.  In  particular,  we  can  take 


Do  =  ^2  S"o  +  ^^4  s; 


(2.2-12) 


which  provides  an  excellent  approximation  of  (2.1-1)  for  any  a 
in  the  aforementioned  range.  The  reason  for  stepping  up  to  the 
fourth  power  in  (2.2-12)  is  to  obviate  any  necessity  of  having 
to  keep  track  of  the  sign  of  So  in  subsequent  applications. 
The  corresponding  expression  for  the  damping  coefficient  C  is 


C 


E 


(  4  ) 


(2.2-13) 


2  2 

where  again  the  condition  C  <  <  I  has  been  applied.  Note  that, 
for  the  case  of  a  =  2  in  the  prior  formulation  (2.2-11),  we 
obtain 


D  -  ttojs  ^2  (2.2-14) 

which  is  the  first  term  of  (2.2-13). 

We  thus  conclude,  that  for  the  one  dimensional  case,  a  modi¬ 
fied  stress  strain  relation  based  on  an  elliptical  hysteresis 
loop,  can  be  readily  constructed  such  that  the  resulting  energy 
dissipation  is  the  same  as  that  given  by  the  empirical  relation 
(2.1-1),  or  its  substitute,  (2.2-12).  This  construction  is  con¬ 
siderably  simplified  by  the  restriction  C*  oif  <<  I  ,  which,  in 
physical  terms,  is  equivalent  to  saying  that*the  stress  amplitude 
is  determined  substantially  by  the  elastic  component  (of  strain) 
only. 

2.3  Biaxial  Stress 

The  method  of  extending  the  stress-strain  relation  (2.2-5) 
to  the  biaxial  case  is  not  obvious.  The  reasons  behind  this  were 
noted  in  the  introduction.  The  most  general  uniaxial  to  biaxial 
stress  state  equivalence  presently  anticipated,  as  expressed  by 
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(1-3),  is  therefore  chosen  for  analysis.  Using  (1-3)  and  (2,1-1), 
the  cyclic  damping  for  biaxial  stress  becomes 


Dq  =  J  (  3  T*  +  A*  I  M  ^ 


while  using  (1-3)  and  (2.2-11),  we  obtain 


D„=  I  J„(  A* 


(2.3-1) 


(2.3-2) 


The  criterion  that  the  stress  amplitudes  remain  sufficiently 
small  so  that  no  effects  of  frequency  or  stress  history  might 
enter  is  now  interpreted  as 


T 


< 


Sl 


(2.3-3) 


Either  of  the  formulations  (2.3-1)  or  (2.3-2)  enable  us  to 
generalize  the  stress-strain  relation  (2.2-5)  to  the  biaxial 
stress  case.  Suppose  we  have  a  plane  stress  problem.  Consider 
the  following  construction  for  the  stress-strain  relations 
written  in  principal  coordinates 


cr 


X 


V€, 


)  +  E  i 


+ 


)  , 


I  -v‘ 


V  €. 


)  +  E*  ( 


(2.3-4) 


The  motivation  behind  the  construction  is  clear:  the  first  terms 
of  (2.3-4)  are  the  usual  relations  of  elasticity  theory,  and  the 
second  terms  are  modeled  on  the  first,  with  constants  E*  and  v  * 
providing  the  analogous  extension  of  constant  C  in  (2.2-5). 

Again,  by  introducing  a  harmonic  straining  action  described  by 


€ 


X 


€ox  COS  (Ug  t 


> 


€y  =  €oy  COStUgf 


(2.3-5) 
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we  can  Integrate  (2.3-4)  for  the  cyclic  energy  dissipation.  The 
result  is 


Do  =  E* 


7rE*Ws(  €q3j  +  2  I/*  €qj^  Cgy  *cy  ^  (2.3-6) 

which  expresses  the  energy  dissipation  in  terms  of  the  ampli¬ 
tudes  of  the  principal  strains. 

We  proceed  next  to  construct  the  alternative  expression  for 
Dg  by  solving  for  the  equivalent  normal  (or  uniaxial)  stress 
and  substituting  this  into  the  empirically  based  expression 
(2.3-1).  Thus  the  stress-strain  relations  may  be  written 

;  (  €ox  "‘oy  ) COSO's*  +  *oy  •  sinwj  t  , 


o-y  =  — t  +E  ojj  (  +-*«„)  sir,  0,^  t  . 

By  subtracting  and  adding  the  two  foregoing  equations,  we  obtain 


<r  = 

y 


-  (  €  —  €  )  cos  to  t  +  U)  (  \  —  V*  )[  €  —  €  )  sin  to  t 

+ ;/  ox  oy  s  s  ox  oy  s  » 


O-  +  (T  =  -  (  €  +  €  )  cos  to  t  +  E’'  to  (  I  +  !/*)(  €  +  €  )  SIO  to  t 

X  y  I  —  1/  ox  oy  s  s  ox  oy  s  » 

from  which  we  can  identify  the  maximum  shear  stress  and  dilata- 
tional  components  to  be 


I,*-  {  (t^) }’‘ox+*oy>  ■ 
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Substituting  these  results  into  (2.3-1)  finally  gives 

+  A*  +(  ]  <  ‘ox-^«cy  >'r. 


(2.3-7) 


The  prior  expression,  given  by  (2.3-6)  is  now  compared  with 
(2.3-7),  from  which  we  obtain  the  condition  that 

^  {  i  [  (  )'  *ox  -  *oy  1“ 

+  A*[(-pr^)  +(E*a,3)“(l  +  r*)“]  (€„+  )']  = 


+  €oy  +€oy  )  . 


(2.3-8) 


A  difficulty  which  now  arises  is  that  the  foregoing  expres¬ 
sion  (2.3-8)  is  the  only  condition  which  emerges  from  this  develop¬ 
ment,  whereas  there  are  now  two,  as  yet  unspecified  damping 
parameters  E*  and  v  in  addition  to  the  special  parameter 
A  *.  Equation  (2.3-8)  therefore  merely  provides  a  constraint 
between  the  three  quantities  E*,  v  *  and  A  *.  In  spite 

of  this  indeterminacy,  several  useful,  specific  formulations  are 
possible  provided  we  accept  certain  restrictive  asstimptions. 

Let  us  consider  the  special  case  where  a  =  2.  Equation 
(2.3-8)  is  then  satisfied  if  we  choose 

xxE*a,3=T{A*[(j§;;)  +(E*u-3ftl+i'*)“]  +  -|- f  +  (  E*(..3ftl-.'*)']  }  , 


E*^-*  =  J  {  A*  [  (  7!^  f  +  ( E*a,3)“(  I  )“]  )%  I-.*  )^]  }_ 


If  we  now  apply  the  restriction  that 


(  I  +  1/ )( I  ~  V*  ) 


(2.3-9) 
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and 


(  I  1/ )( I  +  z/* ) 


<<  I 


the  foregoing  equations  simplify  to  the  following: 


F*= 

^  TTO) 


{  — 

I  l\- 


{\-vf  4(l+z/) 


+  vf  J  ’ 


(2.3-10) 


4A(l+y)^~  3(l~y)^ 
4A*(I+z/)^+  3(1“:/)^ 


It  is  tempting,  at  this  stage,  to  investigate  the  condition 
V  -  V  *,  which  is  strongly  suggested  by  the  way  in  which  the  ori¬ 
ginal  equations  (2.3-4)  were  formulated.  This  condition  gives 


E*  =■ 


A*  = 


3  JE' 


Zv  uj  { \  +  V  )  {  \  —  V  ) 


(2.3-11) 


3  (I  -y  ) 

4  ( I  +y  ) 


where  the  expression  for  A  *  is  particularly  interesting  if  we 
recall  equations  (1-7)  and  (1-10) .  This  indicated  value  of  A  * 
is  1/4  for  the  limiting  case  of  Poisson's  ratio  of  1/2,  which 
checks  with  the  aforementioned  equations,  in  particular,  with  the 
conditions  set  down  for  (1-7).  On  the  other  hand,  for  Poisson's 
ratio  equal  to  1/4,  the  indicated  value  of  A  *  is  3/8,  which, 
although  an  interesting  result,  must  be  treated  as  merely  the  con¬ 
sequence  of  an  expedient  postulate. 

An  interesting  result  may  also  be  obtained  for  the  more 
general  case  where  2  <  O'  <  3,  which  includes  the  entire  range  of 
experimental  values  applicable  to  structural  metals  at  low  stress 
levels.  In  this  case,  in  order  to  avoid  dealing  with  fractional 
values  of  a  ,  we  turn  to  the  expression  given  by  (2.2-12).  This 
two- term  approximation  gives 


D„  =  J-jQ  +  J.Q 


(2.3-12) 


where 


0=  1“ 
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Equating  (2,3-6)  with  (2.3-12),  and  applying  the  restrictions 
given  by  (2.3-9),  we  obtain 

<«ox+‘oy>'} 

{  ^2  +  J’.  [  -f-  ( ifr)  ‘  'ox  -  ‘oy>'  A*  (  (  «<„  +  «oy  >'  ]  } 


(2.3-13) 

Equation  (2.3-13)  is  satisfied,  if,  in  particular,  we  choose 


-L. 

TTCU^ 


{  ^2*+  ^.*14  +2-* 


€  €  “^6 
^ox  ^oy  ^oy 


4  A  (  I  +  1/  )^  -  3  (  I  -  y 
4A^(l  +  y  )^+  3(1-1/ 


(2.3-14) 


where 


(I -I/)  4(1+1/)' 


Once  again,  if  we  look  at  the  special  case  where  v  =  v  * ,  we 
find 


2We(|— l/)(  l+l/)‘ 


,  3  J4 

2(|-i/)(|+i/)’ 


''ox+Z-'ox'oy+'oy’l 


3(l-y) 

^  “  4(l+v)  * 
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It  should  be  noted  that  the  strain  amplitudes  .  ,*ox  and  Cqw 
are  not  necessarily  positive  quantities,  although  the  sign  of  cme 
determines  the  sign  of  the  other  on  account  of  (2.3-5).  However, 
(2.3-6),  which  gives  the  dissipation.  Is  quadratic  In  the  strains, 
so  that  the  sign  Is  unimportant  In  this  case. 


III.  VIBRATING  BEAM  EXAMPLE 


3.1  General  Considerations 

We  consider  the  specific  problem  of  a  cantilever  beam  of  con¬ 
stant  cross  section  which  Is  excited  In  Its  fundamental  flexural 
mode  by  a  forced  periodic  rotation  of  Its  support.  This  problem 
has  already  been  studied  by  Pisarenko  who  used  Davldenkov's  equa¬ 
tions  (1-2)  to  describe  the  beam  material.  An  additional  study  of 
this  same  problem,  but  using  the  stress- strain  relations  proposed 
herein,  is  thus  provided  with  an  interesting  comparison.  For  exam¬ 
ple,  it  can  be  argued  that  the  stress-strain  relations  developed 
in  the  foregoing  section  have  such  complication  that  one  might  as 
well  drop  any  pretense  of  approximate  formulations  and  work  from 
a  much  more  general  development.  It  is  one  of  the  objectives  of 
this  section  to  throw  more  light  on  this  question. 

Pisarenko's  analysis,  which  appears  in  Chapter  Two  of  [2]  , 
uses  asymptotic  expansions  to  generate  successively  higher  ordered 
approximations  of  the  solution  to  the  steady  state  problem.  He 
finds,  however,  that  only  the  first  order  approximation  is  neces¬ 
sary  to  produce  acceptable  accuracy.  It  is  assumed  that  the  same 
situation  holds  true  when  the  new  stress-strain  relations  are  used 
in  place  of  Davldenkov's.  A  comparison  of  results  is  then  possible 
using  first  order  approximations  only. 

3.2  Equations  of  Motion 

The  basic  geometry  and  coordinate  system  for  describing  the 
vibrating,  deformed  beam  is  shown  in  Figure  3,  with  the  following 
notation: 

X  =  horizontal  coordinate  axis  coinciding  with  the 
undeformed,  unrotated  centerline  of  the  beam 

y  =  total  deflection  of  the  centerline  from  the  x-axis 
due  to  rotation  and  deformation 

Q  =  small  angle  of  rotation  of  the  clamped  end  of  the 
beam 

y  =  y  -  Q  Xf  the  deflection  of  the  centerline  due  to 
deformation  from  the  rigidly  rotated  position. 
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This  description  is  based  on  the  concept  of  the  elastic  line  whose 
properties  are  derived  for  a  rectangular  cross  section  of  width  b 
and  height  h  as  shown  in  Figure  4. 

The  appropriate  stress-strain  relation,  as  developed  in 
Section  2,  would  be  based  on  equations  (2.2-4),  (2.2-12)  and 
(2.2-13)  which  give 


(  J-j  +  e'j. 


(3.2-1) 


However,  it  is  instructive  to  base  the  stress-strain  relation  more 
directly  on  the  Davidenkov  description.  Thus  if  we  determine  the 
constants  in  (2.2-12)  by  using  (2.2-3)  where  we  arbitrarily  choose 
m  =  2,  the  corresponding  stress-strain  relation  becomes 


3  TTCUg 


0. 


(3.2-2) 


Using  the  expression  for  infinitesimal  strains  (Figure  3) 

a^y 

€  =  z  — i-  , 
a  X 

(3.2-1)  and  (3.2-2)  become  respectively 


cr  =  E  z 


,  E 

ax"  ^‘*’5 


( 3-2+  E^a-,4  )  3,/^.  }. 


(3.2-3) 


cr  =  E  Z 


dx"  37ru»g  ^0  dt  dx"  ‘ 


(3.2-4) 


The  bending  moments,  M  and  M  for  the  two  cases,  at  any  cross  sec¬ 
tion  are  then  given  by 


4  M 

f  2  2 

-  J  «0 
•^0 


dtdx"  '  '  ’  (3.2-5) 


M  =  El 


a"y  8/xb 


afjLb  r 

TT  •'o 


I  2  y  2  _y 

««  - r  z  dz 

0  °  atax" 


(3.2-6) 


On  neglecting  the  effects  of  shear  deformation  and  rotatory  inertia, 
the  equations  of  motion  are 
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y  =  y  +  0  X  , 

a^y  a^M 

^  at"  ax“  ■  <3-2-7) 

Substituting  (3.2-5)  and  (3.2-6)  into  (3.2-7)  gives  respectively 

^  a  a  X*  ^  a  t  ax^ 


1 

EJi  a  -1 

X  + 

dx*  ^ 

1  + 

Trojg  at  ■’ 

2bJ^ 

d* 

h 

f  * 

a 

TTO^s 

dx* 

Jo 

®  dtdx^ 

z^o'z  =  0  , 


(3.2-8) 


a^  ,  a"0 

at=  "ax^ 


.  8^b 


z^o'z  =  0  , 


(3.2-9) 


as  the  general  equations  of  motion  for  the  elastic  line.  The  cor¬ 
responding  boundary  conditions  are 


X  =  0  , 


X  =  V  . 


(3.2-10) 


The  analyses  of  (3.2-8)  and  (3.2-9)  are  essentially  parallel  to  each 
other:  the  former  (3.2-8)  is  the  one  suggested  for  general  appli¬ 

cation,  while  the  latter  (3.2-9)  is  for  the  specific  purpose  of  mak¬ 
ing  a  comparison  with  Pisarenko's  solution.  The  subsequent  analy¬ 
sis  is  therefore  restricted  to  (3.2-9)  only. 

We  reformulate  (3.2-9)  and  (3.2-10)  by  introducing  the  dimen¬ 
sionless  quantities  x  =  ,  z  =  h  17  ,  and  y  =XVu  ,  where  X  is 

an  arbitrary,  small  parameter.  Furthermore,  on  assuming  that  the 
damping  effect  and  the  excitation,  characterized  in  this  case  by 
the  change  in  the  angle  9  ,  are  both  small,  we  denote 

6  =  6q  cos  cut  , 


(x>y)  =  X  £*  (  ?  ,  u ) . 


The  equation  of  motion  then  becomes 


El 


dC^  K  dtdc' 


7)^  d  Tj  + 


+ 


4 


X  V*  Oq  C  cos  a»  t 


with  the  boundary  conditions 


.(3.2-11) 


u  = 

o 

II 

ot  C  =  0  , 

a"u 

oi 

c 

1 

> 

ot  C  =  1  . 

■ 

»  —  u 

ac’ 

3.3  Method  of  Solution  and  Zeroth  Approximation 


(3.2-12) 


The  method  of  solution  to  be  applied  to  the  foregoing  equa¬ 
tions  is  similar  to  that  proposed  by  Pisarenko.  We  assume  that  the 
excitation  frequency  OJ  ,  is  in  the  neighborhood  of  one  of  the 
linear  natural  frequencies;  say  the  fundamental  natural  frequency 

CL)  ,  and  introduce  the  dimensionless  time 
0 

T  =  cut  +  \//  , 

where  i//  is  an  arbitrary  constant.  Using  the  new  variable  r  , 
(3.2-11)  can  be  rewritten 


^  ^  ^  ac"  Jo  ^0  drac" 

,  =  ^  ^0  (■^)  C  cos  (  T  -  \f/  ),  (3.3-1) 

8ubh®  0 

where  v  =  - ' - :r-  • 

^  vp 

We  now  assume  a  solution  in  the  form  of  the  expansions 

u  =  Uq  +  X  U|  + . 

=  OJq  (  I  +  X  CUq,  + . ) 

(3.3-2) 

0  =  0^  (  I  +  X  0^,  + . )  , 
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where 


u  =  S  (C)  cos  T  . 


(3.3-3) 


On  substituting  (3.3-2)  into  (3.3-1),  we  obtain  the  following 
set  of  linear  differential  equations: 


A  2  ,4  2  ’ 

Ot  p\  Wq  OQ 


(3.3-4) 


El  a^u 


^  =  @0  5  cos  (  T  -  v//  )  -  OJ, 


dr^  plX 


^  dC^  Jo  ~d7d^  ’ 

etc.,  where  cu  =  cug  .  The  corresponding  boundary  conditions  are 


ar,% 

d  u, 


(3.3-5) 


a.  _2, 


d  Ui 


U:  = 


ac 

a^ui  a^ui 


ac^  a^ 


=  0 


=  0 


ot  c  =  0  » 


ot  5=1. 


(3.3-6) 


On  substituting  (3.3-3)  into  (3.3-4)  and  considering  the 
boundary  conditions  (3.3-6),  a  set  of  eigenfunctions 


4>;  = 


J  2  sin  kj  sinh 


—  •[(  cos  kj  +  cosh  kj  )  (  cosh  kj  5  ~  cos  kj  5  )  + 
+  (  sin  kj  -  sinh  kj  )(  sinh  kjC  ~  sin  kj5  )} 


(3.3-7) 

is  generated;  and  the  corresponding  engenvalues  k.  are  detemined 
by  J 


cos  kj  cosh  k j  +  1  =0 


where 


2 

p\  u». 


(3.3-8) 


k|*  = 

J  El 


Hence,  the  zeroth  approximation  is 


2  kjEI 
*^0  =  - 


(3.3-9) 


21 


(3.3-10) 


=  4>JO 


3.4  First  Order  Approximation 
Recalling  that 


h 

V  TTS" 


we  have 


«  =  J —  V  9o  ^ 

^  j."  . 

%  ^  =  ~Y~  "^0  -^  >  0 


(3.4-1) 


If  (3.3-10)  and  (3.4-1)  is  now  placed  in  (3.3-5),  we  obtain 

<^^u,  ,  1  ^  . 

- 7  — 7 - 7“  =  a„  Q  cos  (  r  -  9  )  + 


k;  ac' 


y  h  S 


+  w  S  (h  cos  r  +  ~77r~, —  ^  ko  </)  (^ "  +  <i  "'  )  sin  r 

O'  ^0  32  1  0  ^0  ^0  ^0 


(3.4-2) 


which  is  the  differential  equation  for  the  first  approximation. 

The  corresponding  response  equation  is  obtained  by  multiplying 
(3.4-2)  by  <p  COST  and  <p  sin  r  separately,  and  integrating  each 
result  over  the  length  of  Vne  beam  for  one  cycle.  This  eliminates 
the  time  dependence  and,  of  course,  the  energies  must  balance  over 
the  cycle.  We  thus  obtain 

/o  {  ^0  ^  cos  »//  +  S  =  0 

/I  ^  cs 2 

0  '■  0  ^  32  1  °^o^o  ^0  i  ^0  ^  (3.4-3) 

Carrying  out  the  integrations  and  substituting  for  cUq,  from 
(3.3-2),  finally  leads  to  the  response  equation  for  the  first  order 
approximation . 


er=-T{[ 


4  0o(  cos  ko  +  cosh  ko  ) 
2 

kg  sin  kg  sinh  kg 


r-Ifrl'*'!’  o.‘- 
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where 


T-r 


The  first  order  approximation  of  the  solution  for  the  displace¬ 
ment  is  now  readily  constructed  by  using  the  method  of  Fourier  ex¬ 
pansion.  We  thus  assume  the  solution 


u,  (  ^  j  C  )  =  '/'i  (C)  cos  T  +  (^)  sin  r 

which,  on  subsitution  into  (3.4-2),  gives 

~  =  *<0  {  ^  cos  r//  +  8  0)^,,  } 

ko  {®oC  sinv// 


(3.4-5) 


(3.4-6) 


On  expanding 


v//  =  I  a.  <^. 
'  ] =1  J  J 


=  1  a*4>_ 

'  j=l  J  J 


(3.4-7) 


we  find 


=  4k;eocos^ 

•'o  ^  ^  y 

K  j  Kq  i 


<=- 


/o'  {  ®->  ^  (  Ko  t  '^0 


3.5  Numerical  Result  and  Comparison 

The  result  obtained  in  the  foregoing  section  is  now  applied  to 
the  specific  example  studied  by  Pisarenko.  The  beam  is  made  of 
steel  with  the  following  given  data:  m  =  2,-,  u=  18.6,  b  =  1.12  in. 
h  =  0.59  in.,  \  =  15.95  in.,  E  =  2.96  x  10'  psi.,  and  “IO"^. 

The  plot  of  the  response  curve  given  by  (3.4-5)  is  shown  in 
Figure  5;  also  shown  in  the  same  figure  are  the  corresponding  solu- 
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tlons  up  to  the  first  and  second  order  approximations  by  Pisarenko. 
Obviously  his  second  approximation  blows  up  at  the  peak  point  of 
the  response  curve  where  xj/  =  v/  2  •  This  singularity  can  also 

be  noticed  by  studying  equations  12.4  and  12*5  of 

[3]  ,  which  indicate  that  12,4  cannot  be  satisfied  at\/^— Tr/2 

However,  this  singularity  can  be  removed  if  we  modify  Pisarenko's 
method  by  expanding  the  amplitude  of  the  excitation  instead  of 
the  phase  shift.  This  modification  leaves  the  zeroth  and  the 
first  approximations  unchanged.  Therefore,  although  the  second  ap¬ 
proximation  (by  Pisarenko)  blows  up  at  one  point  it  does  not  imply 
the  invalidity  of  the  first  approximation.  Also,  the  result  far 
away  from  the  singularity,  which  falls  closely  on  that  given  by 
(3,4-4)  might  still  be  used. 

Another  reason  for  the  difference  in  amplitude  and  shape  of 
the  response  curves  between  the  two  cases  is  the  strain  hardening 
effect  of  the  Davidenkov  stress-strain  relation.  This  is  easily 
seen  by  plotting  the  locus  of  the  tips  of  the  hysteresis  loops  as 
shown  in  Figure  6.  The  cyclic  stress  sensitivity  limits  (S.)  for 
some  typical  structural  metals  are  also  indicated  in  the  fi^re, 
beyond  which  values  both  formulations  are  invalid  because  of  stress 
history  effects.  All  of  these  points  are  noted  to  be  below  the 
corresponding  yield  limits,  which  iustifies  the  linear  formulation 
of  the  elastic  component  in  (2.2-5);  otherwise,  as  stated  in 
Section  2,  the  typical  cubic  or  some  other  term  would  have  to  be 
included. 

If  the  difference  in  the  response  (Figure  5)  between  the  two 
cases  is  mainly  due  to  the  method  of  solution  and  magnitude  of  non¬ 
linear  elastic  components,  then  the  use  of  the  elliptical  hysteresis 
loop  should  yield  a  satisfactory  result.  At  least,  the  result 
should  be  as  good  as  that  given  by  Pisarenko's  solution.  This  con¬ 
jecture  is  enchanced  by  plotting  the  two  corresponding  hysteresis 
loops,  given  by  (1-2)  and  (3.2-2)  respectively,  as  shown  in  Figure 
7.  The  strain  amplitude  is  taken  at  a  value  slightly  below  the 
cyclic  stress  sensitivity  limit  of  high  strength  steel.  Neverthe¬ 
less  the  geometric  appearance  of  the  two  loops  appears  to  be  prac¬ 
tically  the  same,  except  for  the  strain  hardening  effect. 


IV.  CONCLUSIONS 


The  analytical  technique  of  describing  material  damping  by  el¬ 
liptical  hysteresis  loops  has  been  applied  to  experimental  damping 
data  to  produce  damped,  stress-strain  relations.  The  objective  of 
this  approach  is  to  produce  an  easily  tractable  mathematical  descrip¬ 
tion  of  the  stress-strain  relation  which  at  the  same  time  adheres 
closely  to  experimentally  observed  behavior.  An  extension  of  this 
development  to  the  biaxial  stress  case  is  also  demonstrated,  although 
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certain  special  assumptions  concerning  the  role  of  dlstortlonal 
and  dllatatlonal  strain  energies  are  necessary  In  order  to  achieve 
this.  An  example  application  to  the  flexural  vibration  of  a 
cantilever  beam  Is  shown  to  Illustrate  the  aforementioned  prop¬ 
erties. 
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